IPMU 12-0233 
ICRR-Report-642-2012-31 



Extended analysis of CMB constraints on 
non-Gaussianity in isocurvature perturbations 



Chiaki Hikage a , Masahiro Kawasaki 6,c , Toyokazu Sekiguchi d and Tomo Takahashi e 

a Kobayashi-Maskawa Institute, Nagoya University, Nagoya 464-8602, Japan 

b Institute for Cosmic Ray Research, University of Tokyo, Kashiwa 277-8582, Japan 

Kavli Institute for the Physics and Mathematics of the Universe, University of Tokyo, 

Kashiwa 277-8568, Japan 

Department of Physics and Astrophysics, Nagoya University, Nagoya 464-8602, Japan 
e Department of Physics, Saga University, Saga 840-8502, Japan 



Abstract 

We study CMB constraints on non-Gaussianity from isocurvature perturbations 
of general types. Specifically, we study CDM/neutrino isocurvature perturbations 
which are uncorrelated or totally correlated with adiabatic ones. Using the data from 
the WMAP 7-year observation at V and W bands, we obtained optimal constraints 
on the nonlinearity parameters of adiabatic and isocurvature perturbations. Our 
result shows that primordial perturbations are consistent with Gaussian ones at 
around 2a level for above mentioned isocurvature modes. 



1 Introduction 



The adiabaticity of primordial density perturbations are one of the most important probes 
to understand the early Universe. Although current cosmological observations suggest 
that primordial fluctuations are almost adiabatic and the contribution from isocurvature 
fluctuations is severely constrained, once any deviation from the adiabatic fluctuations are 
detected even if it is tiny, it readily indicates that there are multiple sources of density 
fluctuations. 

In the view of its importance, there have been many works on constraints on isocur- 
vature fluctuations by using power spectrum [1-10]. However, current cosmological obser- 
vations are so precise that its non-Gaussian nature can also be measured/constrained in 
addition to the power spectrum* 1 . 

Compared to the counterpart of adiabatic fluctuations, observational constraints on 
isocurvature non-Gaussianity using actual data have almost not been studied* 2 , except a 
few works [24, 25] . In Ref . [25] , we have obtained observational constraints on isocurvature 
non-Gaussianity by adopting an optimal bispectrum estimator, focusing on uncorrelated 
type of CDM (baryon) mode. The constraints would give interesting implications to 
particle physics as well since this type of isocurvature fluctuations can be generated in 
some models such as the axion model [26-33] and Affleck-Dine baryogenesis [34-36]. 

In fact, isocurvature fluctuations can also be correlated with adiabatic ones in some 
cases and such correlated type is naturally produced in the curvaton model [37-39] , which 
has been attracting attention due to its possibility of getting large non-Gaussianity. In the 
model, depending on when and how DM or baryon asymmetry are produced, positively 
and negatively correlated isocurvature fluctuations can be generated [19,40-51]. Further- 
more, isocurvature fluctuations can also be associated with neutrino [52] , which generates 
neutrino density isocurvature mode* 3 . Such a neutrino mode can be related to lepton 
asymmetry, thus its constraints may give interesting implications to the lepton sector. 

In the light of these considerations, it would be also worth investigating observational 
constraints on isocurvature non-Gaussianities of correlated type and/or other modes in 
addition to CDM (baryon) one such as the neutrino mode mentioned above. The main 
purpose of this paper is to study the constraints on non-Gaussianity of those kinds of 
isocurvature types/modes by extending the analysis of our previous study [25]. By using 
the optimal estimator which we adopted there, we investigate observational constraints 
on non-linearity parameters for correlated/uncorrelated CDM and neutrino isocurvature 
modes from WMAP7 data. Implications of our results are also discussed for some explicit 

^ 1 Future observations of redshiftcd 21 cm line emission can also probe the isocurvature perturbations 
[11,12]- 

^ 2 See Refs. [13-16] for expected constraints from future observations by using Fisher matrix analysis. 
For theoretical aspects of isocurvature non-Gaussianity, see Refs. [16-23]. 

# 3 There is also another isocurvature mode associated with neutrino called neutrino isocurvature velocity 
(NIV) mode [52], in which there is difference between initial velocity perturbations of neutrino and photon. 
However, there are so far no concrete models which can generate NIV perturbations at super-horizon scales. 
We do not study NIV mode in this paper. 
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model. 



The organization of this paper is as follows. In Section 2, we present the models of 
non-Gaussian isocurvature perturbations which we are to constrain in this paper and their 
CMB bispectra. We also make the Fisher matrix analysis to get some idea for the size of 
the error, which also serves as a check of our results presented in Section 4. Implications 
of our constraints are discussed for the curvaton scenario in Section 5. The final section 
is devoted to conclusion. 

In this paper, we denote the mean value and perturbation of a fluctuating field by bar 
and delta, respectively. A bracket ( ) represents an ensemble average. 

2 Non-Gaussian primordial perturbations and CMB 
bispectrum 

2.1 Non-Gaussian perturbations 

Following mostly Refs. [22,23], here we briefly summarize the formalism to discuss non- 
Gaussianities of primordial curvature and isocurvature perturbations. Based on the 5N- 
formalism [53-56], the curvature perturbation on the uniform density hypersurface of a 
component i, which we denote Q, can be given by 



where N(x,t) is an e-folding number from the initial flat hypersurface to the uniform 
density hypersurface of i at time t. In particular, we denote the curvature perturbation on 
the uniform density hypersurface of the total matter as £ tot . Then we adopt the following 
definition of the curvature ( and isocurvature perturbations Si, 



In this paper, we consider two types of primordial isocurvature perturbations, which are 
CDM one (CI) and neutrino density one (NID). 

In the picture of the separate Universe, the spatial dependence of Ni(x, t) should come 
only from the initial values of the degrees of freedom (or scalar fields) during inflation, 
which we write as <f> a (x) = 4> a + 5(f) a (x) with a subscript a represents a species of scalar 
fields. We can expand the curvature and isocurvature perturbations ( and Si locally in 
terms of 4> up to second order, which leads to 






(2) 
(3) 



C(x, t) = N tot!a (t)S4> a (x) + -N totM (t){5(f) a (x)6Mx) ~ (S<t>a (x) 6<f> b (£))} + • • • , (4) 
Si(x,t) = Si, a (t)5(f) a (x) + -Si >ab (t){5(f) a (x)S(f) b (x) - (5(f) a (x)5(t> b (x))} H , (5) 
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where the coefficients with subscripts a, b, etc. attached to N represent derivatives with 
respect to the initial fields, e.g., 

_ dN tot (x,t) 



(6) 

In this paper, we assume that the initial fluctuations 5<j) a are Gaussian and different 
initial fields are uncorrelated, i.e. (5<f) a 5<j)b) = for a ^ b. Nonlinearity parameters of the 
local type for the curvature and isocurvature perturbations are defined as follows: 

m = cg(x) + pMCo^) - (£(*)», (7) 

S l (x) = S iG (x) + fSl°\sf G (x) - (S? G (x))), (8) 

where ( G and Si G are the Gaussian part of ( and Si, respectively. Eqs. (7) and (8) can be 
recast into a general form: 

X\x) = X£(x) + f^X^xf - (X^(f) 2 », (9) 



where Xq is the Gaussian part of X A . Here we defined a non-linearity parameter 
which is related to the adiabatic and isocurvature ones via f^ L = |/nl and f^ L = f^°^- 

We present constraints in two limiting cases, both of which are of particular interest. In 
one limiting case, Cg an d S iG are assumed to be uncorrelated, i.e. (Cg^g) — 0. This is the 
case that we discussed for CDM (baryon) isocurvature mode in the previous paper [25] . In 
terms of the scalar field fluctuations, this corresponds to the case that the curvature and 
isocurvature perturbations are generated from two different fields 4> and a, i.e. Cg oc 5(p 
and SiQ oc 5a. As we have stated in Introduction, this form of non-Gaussian isocurvature 
perturbations can be generated in e.g. the axion model [26-33] and the Affleck-Dine 
mechanism [34-36]. In another limiting case, £g and Siq are assumed to be totally (anti-) 
correlated, so that (Cg^g) = ±\/ {CgCg) (Si G SiG) . This corresponds to the cases where ( 
and Si G are generated from only a single field 0, so that ( G oc Si G oc 5(f). Such a case can 
be realized in e.g. the curvaton scenario [41,43]. 

From Eq. (9), the primordial bispectrum can be given by 

(X A (h)X B (k 2 )X c (k 3 )) = [f A h P AB (k 2 )P AC (h) + (2 perms)] (2n) 3 5^(k 1 + k 2 + k 3 ), (10) 

where Pab(^) is the power spectrum of (cross-) correlation between Xq and X<f, i.e. 

(X£(k)Xi$)) = P AB (k)(2n) 3 5(k + k'). (11) 

Note that in Eq. (10), we omitted the loop contributions, assuming that the tree terms 
are leading ones in the power spectrum and the non-linearity parameters. 

Current observations have not detected any isocurvature perturbations and only the 
upper bounds on the amplitudes of the isocurvature power spectra have been provided. 
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The upper bounds depend on the kinds of isocurvature modes and their correlations with 
the curvature perturbations. We define the ratio of the amplitude of isocurvature power 
spectrum to that of adiabatic one at the reference scale fc* as 

a= n^T (12) 

In this paper, we take the WMAP convention k* = 0.002Mpc _1 . In the case of CI, WMAP 
7-year result [58] gives upper bounds* 4 

a < 0.15 (uncorrected), (13) 
a < 0.008 (totally correlated). (14) 

For the case of NID, a combined analysis from the WMAP 7- year and the ACT data [59] 
gives upper bounds [60] #5 

a < 0.19 (uncorrected), (15) 
a < 0.013 (totally correlated). (16) 

Here some comments on the parametrization of non-linearity parameters are in order. 
When we consider the bispectrum with multiple sources, one can generally define the 
non-linearity parameters f^ C as follows [14]: 

(X A (k 1 )X B (k 2 )X c (k 3 )) = (2n)W\h + k 2 + k 3 ) 

x ^ L BC P c (fc 2 )P f (A; 3 ) + ^ , 

(17) 

where f^ C = f^ B is satisfied due to the symmetric structure for the latter two indices 
and the (total) curvature power spectrum is used instead of Pab defined in Eq. (11). When 
we consider the case, where curvature and isocurvature perturbations (of a single type) co- 
exist, there exist 6 non-linearity parameters in general: , /£' L S , / N £ C , f^f C , f^f S . 
However, if one considers uncorrelated or totally correlated isocurvature models, these pa- 
rameters can be reduced or related in some way. By using /£ L = §/nl an d /nl = /nl°^ 
defined in Eqs. (7) and (8), for uncorrelated models, the above non-linearity parameters 
are given by 

= /nT = « 2 /nl, f~NL C = (others). (18) 



# 4 In Ref. [58], only the positively-correlated CI model is studied, and does not investigate the negatively- 
correlated one. However, since a allowed observationally is so small, that the upper bounds on a for totally 
correlated cases would not depend much on the sign of the correlation. 

# 5 We note that the effective number of neutrinos iV e ff is also varied in Ref. [60] . However, the degeneracy 
between a and N e a is not severe, thus we expect that constraints with fixed iV c ff = 3 should be almost 
the same. For most recent constraints on a with fixed N c g, we refer to Ref. [1], which uses the WMAP 
3-year data. 
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For totally correlated models, we have 



K,CC _ K K,CS _ ,1/2 K X,SS _ K 

iNL — iNL; JNL — z!za J NL> JNL — U JNL) 

Jnl - "Jnl) Jnl - ±Q; Jnl? Jnl - a Jnl> l iy J 

where the signs should be +/- if isocurvature perturbations are positively/negatively cor- 
related with curvature ones. 

Thus, for uncorrelated modes, 6 parameters reduces to f^ L and the combination 
of a 2 f^ L , which we investigate in the analysis of this paper. For correlated modes, 
we still need 6 parameters. Although in principle, we can probe all of these parame- 
ters/combinations, we focus on f^ L = and «/nl = /nl^ f° r the correlated case since 
the contributions to the bispectrum from other isocurvature modes are weak, hence the 
constraint on an d /nl wou ld be most severe. We will again discuss this issue after 
giving a summary of CMB bispectrum in the next subsection. 



2.2 CMB bispectrum 

Now we briefly summarize the formalism to calculate CMB bispectrum. In the analy- 
sis, we neglect the secondary non-Gaussianities arising from the second or higher order 
cosmological perturbation theory. 

The harmonic coefficients of the CMB temperature anisotropy from primordial pertur- 
bations X A can be given as 



a lm 



/rlk 3 
j^gf(k)x A (k)Y; m (k), (20) 



where gf(k) is the temperature transfer function for the primordial perturbations X A . 
The total CMB anisotropy is the sum of those from the curvature and isocurvature per- 
turbations, i.e. ai m = a< !m + a ? m - The angular power spectrum Ci, defined by (ai m ai' m *) = 
Ci^ii'^mm', can be written as 

Q = -J2 [ dkk 2 gf(k)gf(k)P AB (k). (21) 

^ AB 

It is convenient to define the reduced bispectrum b^i^ as 

(ai (22) 

where ^ 3 2m3 is the Gaunt integral, which can be written in terms of the Wigner-3j 
symbol as 



G i lhh _ / (2/ 1 + l)(2/ 2 + l) (2/ 3 + l) / h l 2 l 3 \ (h k h , 
y mim2m3 y 4?r I mi m 2 m 3 I 1 0'' 
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Given the bispectrum of primordial perturbations in Eq. (10), the reduced CMB bispec- 
trum is 

b hhh = 2 /nl / drr 2 (a£(r)^(r)^(r) + (2 perms)) , (24) 



ABC 



where af(r) and /3 t ' (r) are 



aftr) = -J dkk'gfWUkr), (25) 

/f"V) = IJ dktfgfWJiWPABik). (26) 

By defining 

bfff 3 = 2 1 rfrr 2 (<(r)/^ A (r)/3^(r) + (2 perms)) , (27) 
Eq. (24) can be recast into 

&ii«s = £[/&^ + (2pennfl)]. (28) 

ABC 

For later convenience, we introduce a normalized bispectrum 

?A _ 9b h i 2h 
U J] 



NL 



^ZKhh + ( 2 Perms)], 



with which we can rewrite Eq. (28) as 

&IiW» = £/m#W.- (30) 

When we consider primordial perturbations of uncorrelated isocurvature type, we ob- 
tain 



f ^S,SS 
1/2^3 hhh' 



Since &^/ 2 j 3 oc a 2 /^ L , CMB anisotropy is affected by the combination of a 2 f^ L . Thus what 
we can constrain from CMB experiments is ck 2 /n L , not alone. 
For the case of correlated isocurvature perturbations, we obtain 

tfihh = ^hhh + ^hhh + ^hhh-> ( 33 ) 

is _ tS.CC , 2h S ' < > S 4- h S ' SS (M) 

°hhl 3 - °hhh + Z °hhh + ^iWs" ^ 
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In terms of a, the second and third terms in the right hand side in each equation are 
respectively of 0(\fa) and 0(a) higher order than the first terms. From current observa- 
tions, a should be small as a < 0.01, hence the first terms in Eqs. (33) and (34) are leading 
in the total bispectrum of Eq. (30), which are proportional to /^ L and a/ff L , respectively. 
Thus the constraints on these quantities are most severe among the non-linearity param- 
eters/combinations listed in Eq. (19). However, as mentioned in the previous subsection, 
observations are in principle also sensitive to the amplitudes of &ff 2 j 3 , %z 2 z 3 j and 
bf^ij , which are respectively proportional to a 1//2 /^ L , af^ h , a 3 ^ 2 f^ L and « 2 /nl- Although 
their contribution could also affect constraints on / NL and af^°\ nevertheless, such an 
effect is expected to be weak as far as a is small as a < 0.01. We will come back to this 
issue at the end of Section 4. 



2.3 Fisher matrix analysis 

Before presenting our constraints from observed data, we here present expected errors 
from WMAP7 data based on the Fisher matrix analysis, which serves as a cross-check of 
our results given in Section 4. 

According to Refs. [61,62], given the bispectrum of Eq. (30), the Fisher matrix for the 
bispectrum is given by 

rp /sky V- (2h + l)(2l 2 + 1)(2Z 3 + 1) (h hh\ 2 KlilMhh (V * 

~q 2. ^ [ooo c h c h c l3 ' ( } 

where Q = Ci + Ni with C\ and Ni being the power spectra of signal and noise, respectively. 
Given a survey with a finite resolution, Fab saturates as / max are increased. 

The signal power spectrum C\ are computed using the CAMB code [63] assuming 
the mean parameters from the WMAP 7-year data alone [58]. We set the same spectral 
indices for both the power spectra of adiabatic and isocurvature perturbations. On the 
other hand, the noise power spectrum can be approximated with the Knox's formula [64] 
as, 



N i = ^fwhm^t ^xp 



(36) 



81n2 

where #fwhm is the full width at half maximum of the Gaussian beam, and ot is the root 
mean square of the instrumental noise par pixel. The total Ni of a multi-frequency survey 
can be given by a quadrature sum of Ni of each frequency band. For the Fisher matrix 
analysis, we adopt the specification of the 7-year WMAP V and W bands for 6*fwhm and 
<tt, which are are listed in Table 1 [65]. Furthermore, we assume that / s k y = 0.72, which 
is consistent with the mask we adopt in Section 4. 

From Fab, we can evaluate ideal constraints on the nonlinearity parameters. For a 
single nonlinearity parameter, two different constraints can be derived. One is the non- 
marginalized constraint in which we assume that the other nonlinearity parameter is fixed 
to zero. In this case, 1 a error should be given by A/^ L = I/^/Faa- Another one is the 
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band 


V 


W 


#fwhm [arcmin] 


21.0 


13.2 


o T [mK] 


25 


38 



Table 1: Survey parameters for 7- year observation of the WMAP V and W bands. 



marginalized constraint which assumes the other nonlinearity parameter varies freely. In 
this case, la error is given by A/^ L = J F^\- In Table 2, we list ideal la constraints 
from the WMAP 7-year observation for uncorrelated CI and NID models. Those for the 
correlated CI and NID models are listed in Table 3. These results are to be consulted for 
validity check of our constraints from actual data in Section 4. 



uncorrelated CI 20 (22) 65 (72) 
uncorrelated NID 20 (38) 135 (256) 

Table 2: la errors evaluated from the Fisher matrix analysis on the nonlinearity pa- 
rameters from the WMAP 7-year observation for the uncorrelated CI and NID models. 
Errors with and without parenthesis correspond to non-marginalized and marginalized 
constraints, respectively. 



correlated CI 20 (24) 120 (140) 
correlated NID 20 (56) 51 (143) 



Table 3: The same as in Table 2 except for the correlated CI and NID models. 



3 Analysis method 

Our analysis method is the same as in our previous paper [25]. For more details, we refer 
to Ref. [25] and references therein. 

We adopt the following optimal estimator of nonlinearity parameters f^ L , 



Ml\ _ | ^P rim ^/NL =1 ^Prim)^' 
Ml/ \(^prim)/£ =1 




(37) 
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where a square bracket with subscript f$ L = 1 indicates an ensemble average over simu- 
lations with non-zero non-linearity parameter shown in the subscript and the other non- 
linearity parameter being fixed to zero. Sp rim is the cubic statistics [66-69] , which is given 

*Sprim = g ^ ] ^ 1 i 2 «3^mi 2 m 2 m3 [^hmi ^Z 2 m 2 ^3^13 ~ 3a; imi (5; 2m2 a^ 3m3 )o] , (38) 
{lm} 

where the square bracket with subscript indicates an ensemble average over Gaussian 
simulations. a/ m is a harmonic coefficient of observed (or simulated) data maps weighted 
by inverse- variance, 

aim = y^^imMm'^'m', (39) 
I'm' 

which is derived from the optimality of the estimator of Eq. (37) [70]. Eq. (37) can be 

schematically represented as f^ L = ^^(£p r i m )/B =1 Sp Tim - In our analysis, to evaluate 

b NL 

(aimO'i'm')o, we accumulated at least 250 Monte Carlo (MC) samples. 

'NL? 



Assuming the Gaussianity for X A , we estimate the covariance of the estimator f A 



(/nl/nl)o ^ ] (^prim) fA' =1 (^prim^prim)o(^prim) f& =l 

(40) 



A'B' 



NL JNL" 



\ prim/ jB = -^- 

Here we used the relation (S A lixa ) fB =1 = (S A im Sp Tixa )o, which has been proved in our 
previous paper [25] . When we assume that either of adiabatic or uncorrelated isocurvature 
perturbations are non-Gaussian and the other is Gaussian, Eq. (37) can be reduced to 

/nL = ^prim/ {^piim) / N l= 1 ' (^) 

where the subscripts A indicates the perturbations which are assumed to be non-Gaussian. 
In this case, the variance of /^l i s gi ven by 1/ (5^. im ) =i- 

In order to determine the normalization (S^^fB _•, , we need to simulate non-Gaussian 

\ prim/ J^L - 1 

CMB maps. We adopt the method for local-type non-Gaussianity developed in Ref. [71], 
in which the quadrature method with optimized nodes and weights for line of sight integral 
is used. In our analysis, we require the method to be accurate such that the mean square 
of the error is less than 0.01 at each multipole (lm). For Z max = 1024, we found that this 
level of accuracy requires 42, 15 and 14 nodes for adiabatic, CI and NID perturbations, 
respectively. 

Computation of inverse-variance weighted map a\ m is performed based on the method 
of Ref. [72], which adopts a conjugate gradient method with multi-grid preconditioning. 
This method also allows us to marginalize over the amplitude of components which are 
proportional to template maps. 
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4 Result 



In our analysis, we assume a flat power-law ACDM model as a fiducial one and adopt the 
mean values for cosmological parameters from the WMAP 7-year data alone [58] 

(tt b , n e , H , t, A adi ) = (0.0448, 0.220, 71, 0.088, 2.43 x 10~ 9 ). (42) 

Here and Q c are the density parameters for baryon and CDM, respectively. Hq is 
the Hubble constant in units of km/sec/Mpc and r is the optical depth for reionization. 
A a di is the amplitude of power spectrum of adiabatic perturbations at the reference scale 
K = 0.002MPC- 1 . With this definition, the power spectrum can be written as V^(k) = 
(k 3 /27i 2 )P^(k) = Aadi(A;/A;*) nadi_1 , where n adi is the spectral index for P^. We assume 
that the power spectrum of isocurvature perturbations is also given in a power-law form 
"PsgSg(^) = (^ 3 /2 7I " 2 )-Ps , g 5 , g(^) k ^™ iso ~ 1 - With regard to the values of n a di and n iso , 
we adopt several different settings, which will be separately presented when we show 
our constraints. The transfer function of CMB is computed using the CAMB code [63]. 
We combine the foreground-cleaned maps of V and W bands of the WMAP 7-year data 
[73, 74]* 6 with a resolution N B ^ e = 512 in the HEALPix pixelization scheme [75]* 7 . We 
adopt the KQ75y7 mask [74], which cuts 28.4 % of the sky. We also set the maximum 
multipole Z max to 1024 in our analysis. We marginalize the amplitudes of the monopole 
I = and dipoles I = 1 as default and also optionally marginalize the amplitudes of the 
Galactic foreground components at large angular scales using the synchrotron, free-free 
and dust emission templates from Ref. [74]. 

Constraints for the case of uncorrelated isocurvature perturbations are listed in Table 4. 
Concerning the spectral indices, we fix n a di to the WMAP 7-year mean value 0.963, while 
for n iso , we adopt two settings, n iso = 0.963 and n iso = 1, which are separately analyzed. 
Joint constraints on / NL = §/nl anc ^ a2 f^i°^ = ^Vnl are shown in Fig. 1. First we found 
that the la errors are roughly the same as those obtained from the Fisher matrix analysis 
given in Table 2 regardless of marginalization of foregrounds and fiducial nj S0 . The fiducial 
values of nj SO little affect the central values for both the uncorrelated CI and NID models. 
While the marginalization of foregrounds can bias the central values in the uncorrelated 
CI model by as much as 0.5<r, the effects in the uncorrelated NID model are smaller, from 
which we can see that the effects of residual foregrounds are moderate. To conclude the 
case for the uncorrelated modes, the primordial perturbations are consistent with Gaussian 
ones at around 2a level, even if we consider uncorrelated CI or NID perturbations. 

Constraints for the case of correlated isocurvature perturbations are listed in Table 4. 
With regard to the values of spectral indices, we adopt two different cases, n adi = n iso = 
0.963 and n adi = n iso = 1 for correlated case. Joint constraints on / NL = §/nl an< ^ 
a f^i°^ = a Mh are shown in Fig. 2. We found that the resulting constraints are consistent 
with the Fisher matrix analysis shown in Table 3 in this case as well, although marginalized 

# 6 http://lambda.gsfc. nasa.gov 
# 7 http: / /healpix.jpl. nasa.gov 
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setups 


/nl 


a 2 f (iSO) 
« Jnl 


CT n- — 




7 / r\ i~ m t~\ '1 1 - o v > i q t*o*i n q 1 i7Q linn 
W/U LClilUlcLLt; lilcxl^lllcLllZjcLljlUll 


(50 ± 23) 


±0 nz uu 

(-51 ±72) 






WT tpmTllafp TCnn 3 1 T7£lf 1 on 
W / udllJJlCXljC lllCXl bllllCXllZiClulWll 


37 + 91 
(41 ± 23) 


99 + 64 
(-28 ±71) 




1 


W/U LtslilUlcLL" lIlcLl^lllcLllZjdLlUll 


(51 ±23) 


(-34 ±69) 






w / fpmTilafp ma ycnn a 1 T7fif 1 on 

W / udllJJlCXljC lllCXl bllllCXllZiClulWll 


33 ± 91 
(35 ± 23) 


30 ± 66 
(-15 ±72) 


NID, n iso 


= 0.963 


w/o template marginalization 


43 ± 21 
(65 ± 39) 


191 ± 140 
(-173 ±261) 






w/ template marginalization 


34 ±21 
(48 ± 39) 


164 ± 143 
(-116 ±266) 


NID, n iso 


= 1 


w/o template marginalization 


40 ±21 
(57 ±40) 


178 ± 137 
(-133 ±257) 






w/ template marginalization 


36 ±21 
(48 ± 40) 


175 ± 137 
(-87 ±257) 



Table 4: Constraints on / NL and a 2 / NL at la level for the cases of uncorrelated isocur- 
vature perturbations. A value with (without) parenthesis is a constraint on a nonlinearity 
parameter without (with) marginalization of the other one. 



errors for the NID model are substantially larger than the Fisher analysis, which may 
result from the inhomogeneous noise and the sky cut. We found that fiducial values of the 
spectral indices little affect the constraints. On the other hand, template marginalization 
can bias the central values by as much as 0.5cr, in particular for CI model. However, as in 
the uncorrelated case, the effects of residual foregrounds are regarded not to be severe at 
all. We can conclude that the primordial perturbations are consistent with Gaussian ones 
at around 2a level also in the case of correlated CI and NID models. 

Finally we comment on the contributions from the terms which are higher order with 
(D(y/a) and 0(a) such as the second and third terms in Eqs. (33) and (34). In the analysis 
of correlated isocurvature models, we have neglected those terms since they are expected 
to give little effects on the constraints. To check this, we have also made an analysis where 
we included those terms. We have compared the constraints with and without those terms 
for the case of a = 0.003. It is found that the errors are almost unchanged and the central 
values change by no more than O(0. l)a, which justifies ignoring higher order terms in a 
in our analysis. 
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Figure 1: Joint constraints on /nl and ck 2 /nl a ^ 1 anc ^ 2cr levels for uncorrected 
isocurvature models. Different panels show cases for different isocurvature modes and 
values of n iso . Left and right panels respectively correspond to the cases of CI and NID 
modes. Top and bottom panels respectively correspond to the cases of n iso = 0.963 and 
n iso = 1. Here n a di is fixed to 0.963. In each panel, solid red (green dashed) lines are the 
constraints with (without) template marginalization of foregrounds. Blue star corresponds 
to the case where primordial perturbations are purely Gaussian. 
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Figure 2: Joint constraints on /nl and af^ L at 1 and 2a levels for correlated isocur- 
vature models. Left and right panels show the constraints for the cases of CI and NID 
modes. We here assume n a di = n; so and cases of n a di = n iso = 0.963 and 1 are respectively 
shown in top and bottom panels. Other conventions are the same as in Fig. 1. 
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n Qfi^ 
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W/U LtJlllUlClLtJ lilcLl gjlllcLllzjcLLlUll 


41 4- 91 
(50 ±25) 


7fi 4- 1 1 4 

(-82 ± 138) 






W/ tPTYlTHPfp TTIPTCflTl P ll y^flOTl 
W / UClll UlCX LC- lllCXl tilllCXllZiClulWll 


34 ± 91 
(37 ±25) 


90 ± 1 90 
(-26 ± 144) 


( 1 />n . '/-■) 

^J-) "-iso "-adi 


1 


~\~K7 / r\ f~ o tyi r~\ -i I - <:_> 1 1 1 o t* m n a i7Qf i /~\n 
W/U LClilUldLt; lilcLl ^lllcLllzjdLlUll 


4D 4- 91 
(48 ± 25) 


7H _(_ 1 1 4 
(-79 ± 138) 






W/ fPTYlTllafp TTI^TCflTl 3 ll 1 OH 
W / UdllJJlCXLC' lllCXl cllllCXllZiCXljlWll 


37 ± 91 
(40 ± 25) 


qq 4- 1 1 7 
(-25 ± 141) 


NID, riiar, = Tl a A\ 
5 "ISO cAQl 


= 0.963 


w/o template marginalization 


45 ± 21 
(93 ± 86) 


103 ± 55 
(-126 ±220) 






w/ template marginalization 


35 ±21 
(55 ±80) 


82 ±54 
(-53 ± 203) 


NID, n iso = n adi 


= 1 


w/o template marginalization 


42 ±21 
(72 ± 75) 


99 ±53 
(-78 ±191) 






w/ template marginalization 


36 ±21 
(67 ±80) 


86 ±53 
(-80 ± 204) 



Table 5: Constraints on / NL and af^ L for the cases of correlated isocurvature pertur- 
bations. 



5 Implications to the curvaton scenario 

In this section, we demonstrate an application of our result to a model of the very early 
Universe based on the curvaton scenario. We here focus on a model where inflaton and 
curvaton fields source the primordial curvature and CDM isocurvature perturbations. In 
addition, we assume that some fraction of CDM is generated before curvaton decay, and 
the rest of CDM is generated directly from the curvaton decay. No CDM production 
after curvaton decay is assumed. Here we briefly outline the model and its predictions 
for the power spectrum and non-linearity parameters. For more details, we refer to e.g. 
Refs. [18,20-22]. 

Adopting the sudden decay approximation, the curvaton decays synchronously on the 
uniform-density hypersurface of the total matter. On this hypersurface, the energy con- 
servation of the total matter and CDM should be given by 

p r = pWe 4 « W -V+p a e^-V, (43) 
PCDueXfox-O = piV 3(CW) " C) + ptcV 3 ^- , (44) 

where pi^ and p a are respectively energy densities of radiation and curvaton just before 
decay. The subscript <\> for py 1 is put to emphasize that the radiation originates from the 
inflaton decay (or decay products). p T and Pcdm are respectively ones for radiation and 
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CDM after the decay. p^M an< ^ Pcdm are respectively energy densities of CDM produced 
from the inflaton and the curvaton, which leads to Pcdm = Pcdm + Pcdm- * s the 
curvature perturbation generated from the inflaton. ( a and Ccdm respectively are the ones 
on uniform curvaton and CDM energy density hypersurfaces. ( is for the total one. 

In the first equation, we assumed that CDM is energetically negligible during the decay 
and also the Universe is almost completely dominated by radiation after the curvaton 
decay. Otherwise, the Universe is dominated by matter soon after its decay, which would 
require further entropy production. 

To describe the model, we introduce some parameters. We define the fractions of 
energy densities created from the curvaton decay in radiation and CDM as e r and €cdm, 
respectively, i.e. 

Cr = Pa/Pr, e CDM = PcDm/PcDM- (45) 

Then Eqs. (43) and (44) can be recast into 



'1 - e r )e 



4(C(«-C) 



+ e r e 



3(C CT -C) 



3 3(CcDM-C) 



;i-6cDM)e 3 ^ ) -« + eC DMe 3 ^-«. 



(46) 
(47) 



By solving Eqs. (46) and (47) for the curvature perturbation ( and CDM isocurvature 
fluctuations Scdm = 3 (Ccdm — C) U P to second order in £W an d S a = 3(( a — C^)j we 
obtain 



c = c w + 



rS a (1 — r)(3 + r) / rS a 



+ 



2r 



5, 



CDM 



Ccdm 



r 



rS n 



+ 



ecDM(l-ecDM) (l-r)(3 + r) 



2r 2 



6r 



where r is defined as 



3e T 



(rs a y 



(4? 



(49) 



(50) 



Furthermore, can be expanded up to the second order as, assuming that the curvaton 
has a quadratic potential, 

flt = 2^ -(**]", (51) 



a* 



a* 



where a* = a* + 5a* is the value of the curvaton at the time of horizon exit and its 
fluctuations 5a* is assumed to be Gaussian. We also treat £w as a Gaussian field. 

Although our setup here a bit differs from one in Section 2, where curvature and 
isocurvature perturbations are assumed to be generated from a single source, our result 
in Section 4 for correlated isocurvature models can be still applicable as bispectra are 
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generated from the single source S a . Then a and the parameters characterizing amplitudes 
of bispectra, /nl and afj^°' , are given by 



a. 



ccdm 



NL 



9 
A 

5(1 -r)(3 + r) 



r 



.(ISO) 



9 



Q rA 2 

ecDM(l-ecDM) (l-r)(3 + r) 



2r 2 



6r 



(52) 
(53) 
(54) 



where we defined a parameter A, which is the ratio of total curvature power to the con- 
tribution of the curvaton therein, i.e. 



A = 



(C 2 ) 



<(rS„/3) a >- 



(55) 



Fig. 3 shows constraints on r and ccdm at 2a level with several fixed values of A. 
Shaded regions are allowed by the constraints a, /nl and atf^°\ Regarding the upper 
bound on a, when A>1, Scdm is almost uncorrelated with the curvature perturbation 
since it almost originates to that from the inflaton. Thus Eq. (13) is adopted. As for 
bounds on /nl and f$i°\ we adopt those for correlated models with n at a 
template marginalization, /nl < 90 and —307 < otfm, < 257. 



1 and 



When A is as large as (9(100), there is some parameter region where ck/nl can be 



observably large as ck/nl^ 



case, the constraints on af^°^ becomes important 



0(100) without producing too large a and /nl- In such a 



6 Conclusion 

In this paper, we have investigated non-Gaussianties in isocurvature fluctuations by ex- 
tending our previous analysis for CDM uncorrelated mode [25] to a more general type. 
We derived optimal constraints on them using the CMB temperature maps of the WMAP 
7-year observation. Specifically, we have studied correlated/uncorrelated CDM and neu- 
trino density isocurvature modes along with the counterpart for the adiabatic one. Their 
constraints are summarized in Tables 4 and 5 for uncorrelated and correlated modes, 
respectively. We have also presented joint analysis in which both the adiabatic and isocur- 
vature fluctuations are assumed to exist and obtained 2D constraints as shown in Figs. 1 
and 2 for uncorrelated and correlated cases, respectively. We have checked our results by 
making the Fisher matrix analysis adopting the specifications of WMAP7 experiments, 
which should give an idea how large the errors can be. Comparing with the results ob- 
tained by using actual CMB temperature maps, we have found that they are almost the 
same in size, which shows the reasonability of our constraints. 

Our results indicate that both of the adiabatic and isocurvature perturbations, even if 
we consider various types for the latter, are consistent with Gaussian ones at around 2a 
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Figure 3: Constraints on curvaton scenario in r-ecDM plane with fixed values of A. From 
top left to bottom right panels, A scales from 10 to 10 4 . An arrow indicates the direction of 
which region is allowed for each constraints. Shaded regions are allowed by the constraints 
oi, /nl and af^°\ For the case of A = 10, the region allowed by the constraints from a 
is the one between two red lines. See text for more details. 
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level. Although we have not detected any non-Gaussian signatures in those perturbations, 
their constraints may give interesting implications to some explicit models. As such an 
example, we have considered a curvaton scenario where correlated CDM fluctuations are 
generated. As mentioned in the text, isocurvature fluctuations have already been severely 
constrained by power spectrum measurements. However, we argued that, in some par- 
ticular cases, our constraints for isocurvature non-Gaussianities can give a more stringent 
one. 

In the near future, we are expecting to obtain much more accurate measurements of 
CMB from Planck. At this era of precision cosmology, isocurvature fluctuations and their 
non-Gaussianities would provide fruitful information on the physics of the early Universe 
and the research along this line may also give important implications to particle physics 
as well. 
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